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摘 要 : 为 了 提高 过 感 图 像 的 实时 分 类 准确 率 与 效 窑 ， 提 出 了 一 种 基于 蚁 群 优化 算法 与 独立 特征 集 的 阔 感 图 像 集 实 


时 分 类 算法 。 首 先 ， 提 取 帝 感 图 像 的 小 波 域 特征 与 颜色 特征 ， 并 且 组 成 特征 向 量 ; 然后 ， 采 用 蚁 群 优化 算法 对 特征 
空间 进行 优化 ， 独 立地 选 出 每 个 分 类 的 显著 特征 集 ， 从 而 降低 每 个 子 特 征 空 间 的 维度 ; 最 终 ， 每 个 分 类 独立 地 训练 


一 个 极限 学 习 机 分 类 器 ， 从 而 实现 对 遥感 图 像 集 的 分 类 。 基 于 公开 的 须 感 图 像 数据 集 进 行 了 仿真 实验 ， 结 果 显示 本 
算法 实现 了 较 高 的 分 关 准 确 率 ， 并 且 实 现 了 较 高 的 计算 效率 。 
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Abstract: In order to improve the accuracy and efficiency of real time classification of remote sensing images, a real time 
classification algorithm of remote sensing images based on the ant colony optimization algorithm and independent feature 
sets is proposed. Firstly, wavelet features and color features of remote sensing images are Abstract: d, and the features form 
the feature vectors; then, the ant colony optimization is adopted to optimize the feature space, and the significant feature set 
of each class are selected independently to reduce the dimension of each feature sub-space; lastly, an independent extreme 
learning machine is trained for each class to realize the remote sensing images classification. Simulation experimental 
results based on the public remote sensing image dataset show that the proposed algorithm realizes a good classification 
accuracy and computational efficiency. 

Key words: artificial intelligence; feature abstraction; remote sensing image; computational efficiency; ant colony 
optimization algorithm; extreame learning machine 


0 ”引言 然后 将 这 两 个 提取 的 特征 串 行 拼接 ， 该 方案 使 用 使 用 支持 向 
于 量 机 作为 遥感 图 像 集 的 分 类 器 。 该 方案 有 效 地 提高 了 图 像 分 

遥感 技术 是 军事 侦察 、 导 弹 预警 、 军 事 测 绘 、 海 洋 监视 类 准确 率 ， 但 是 时 间 效 率 较 低 。 文 献 [9] 提 出 了 一 种 基于 显著 
与 气象 观测 等 领域 的 重要 研究 课题 之 一 ， 其 中 通过 遥感 图 像 ” 特征 和 GVF Snake 的 高 分 辩 率 遥感 图 像 道路 提取 方法 ， 通 过 
对 地 面 进行 探测 与 侦查 具有 巨大 的 研究 价值 由 。 在 矿产 检测 融合 颜色 对 比 度 和 空间 统计 特征 计算 显著 性 图 ， 并 以 输出 的 
与 石油 勘探 工程 应 用 中 ， 通 过 卫星 遥感 图 像 对 地 表 进 行 初步 显著 图 最 大 值 作为 GVF Snake 的 初始 种 子 点 。 该 方案 同时 提 
hs 高 了 计算 效率 与 检测 精度 ， 但 是 由 于 该 方案 采用 了 GVF 
在 机 器 人 探索 芭 漠 地 区 的 工程 应 用 中 ， 采 集 遥 感 图 像 也 是 探 snake 的 种 子 生 长 算法 ， 导 致 该 方案 仅 对 轮廓 显著 的 目标 具 
测 机 器 人 的 一 个 重要 任务 B1。 但 是 探测 机 器 人 与 卫星 的 通信 有 效果 。 
能 力 有 限 ， 如 果 将 遥感 图 像 集 实时 地 发 送 至 工作 站 ， 则 会 站 小 波 特征 能 够 提取 遥感 图 像 的 细节 信息 00， 本 文 则 同时 
来 巨大 的 通信 负担 内。 因此 ， 在 探测 机 器 人 或 者 卫星 中 实时 考虑 了 遥感 图 像 的 小 波 域 特 征 与 颜色 特征 。 为 了 保证 较 低 的 
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地 对 般 感 图 像 集 进行 预 处 理 ， 不 仪 可 以 减轻 通信 负担， 而 遥感 图 像 分 类 处 理 时 间 ， 有 效 地 减少 特征 数量 则 是 另 一 个 重 
出 有 助 于 对 地 表 的 定位 效果 59。 要 的 课题 。 目 前 ， 卓 有 成 效 的 特征 优化 方案 主要 有 蚁 群 优 化 

遥感 图 像 的 细节 信息 丰富 、 特 征 维度 高 ， 一 般 采 用 传统 ”算法 00、 粒 子 群 优化 算法 02、 遗 传 算法 03 等 。 其 中 蚁 群 优化 
的 特征 提取 方案 提取 遥感 图 像 的 特征 ， 获 得 的 特征 表达 能 算法 取得 了 更 好 的 效果 ， 因 此 本 文 也 采用 蚁 群 优 化 算法 对 每 


弱 并 且 信息 损失 严重 ， 对 物 感 图 像 后 续 的 处 理 具有 不 利 的 影 ”个 分 类 的 特征 集 进行 优化 处 理 。 分 类 器 的 时 间 效 率 也 是 决定 
响 中 。 文 献 [8] 对 预 处 理 后 的 图 像 进行 降 维 处 理 ， 提 取 像素 的 。 ”遥感 图 像 分 类 效率 的 一 个 重要 因素 ， 目 前 遥感 图 像 分 类 领域 
光谱 维特 征 ,并且 使 用 卷 积 神经 网 络 提取 图 像 块 的 空间 特征 ， 。 中 常用 的 分 类 器 主要 为 支持 向 量 机 n4， 本 文通 过 实验 发 现 极 
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限 向 量 机 的 计算 效率 略 优 于 支持 向 量 机 ， 因 此 本 文采 用 极限 
向 量 机 作为 分 类 器 。 

综 上 所 述 ， 为 了 在 遥感 测量 终端 实现 对 遥感 图 像 集 的 实 
时 分 类 处 理 ， 应 当 在 保证 遥感 图 像 分 类 准确 率 的 情况 下 ， 有 
效 地 降低 算法 的 计算 时 间 。 本 文 提 出 了 一 种 基于 蚁 群 优化 算 
法 与 独立 特征 集 的 遥感 图 像 集 实 时 分 类 算法 ， 该 算法 提取 遥 
感 图 像 的 小 波 特 征 与 颜色 特征 ， 采 用 蚁 群 优化 算法 对 特征 数 
量 进行 优化 处 理 ， 为 每 个 图 像 分 类 分 别 训练 一 个 独立 的 极限 
学 习 机 分 类 器 。 最 终 基 于 公开 的 遥感 图 像 数据 集 进 行 了 仿真 
实验 ， 结 果 表 明 本 算法 在 保证 较 高 图 像 分 类 准确 率 的 同时 ， 
有 效 地 提高 了 分 类 处 理 的 计算 效率 ， 有 助 于 在 遥感 测量 终端 
实现 实时 的 图 像 预 处 理 。 


1 ”特征 提取 


寺 征 提取 从 复杂 的 图 像 数 据 集 中 提取 出 一 个 特征 向 量 ， 
采用 特征 向 量 表征 一 个 大 规模 的 数据 集 ， 由 此 节约 数据 的 存 
储 成 本 与 计算 成 本 。 在 本 文 的 特征 提取 实验 中 ， 从 每 个 像素 
的 窗口 提取 该 像素 的 特征 向 量 。 地 表 遥 感 图 像 中 的 岩石 特征 
采用 小 波 变换 中 的 纹理 特征 描述 符 表示 ， 将 小 波 域 的 统计 参 
数 作为 像素 的 描述 符 。 其 他 的 地 表 特 征 则 采用 颜色 特征 描述 
符 表示 。 本 文 的 特征 描述 符 共 包括 四 种 特征 ， 分 别 为 : 小 波 
域 的 系数 特征 、 主 颜色 描述 符 (DCD) 特 征 、 局 部 颜色 直方 图 
特征 与 颜色 统计 特征 。 
1.1 小 波 特 征 
首先 ， 为 输入 图 像 的 每 个 像素 建立 一 个 窗口 ， 对 窗口 进 
行 一 级 的 二 维 小 波 分 解 处 理 。 然 后 ,将 欧式 范 数 分 别 放 入 LH、 
HL 元 素 各 行 、 各 列 的 小 波 系数 中 05。 最 终 ， 分 别 计 算 欧 式 
范 数 的 平均 值 与 方差 ， 获 得 8 个 小 波 特征 值 。 
1.2 颜色 特征 

颜色 特征 属于 低级 别 特征 ， 颜 色 特 征 对 于 旋转 、 缩 放 等 
变化 不 敏感 ， 遥 感 图 像 存 在 较 多 的 噪声 ， 并 且 许 多 图 像 存在 
旋转 、 缩 放 等 变形 ， 因 此 颜色 特征 可 以 解决 这 些 问题 。 本 文 
采用 了 三 个 颜色 特征 , 分 别 为 : 主 上 颜色 描 述 符 (dominant color 
descriptor,DCD) 特 征 、 局 部 颜色 直方 图 与 颜色 统计 特征 。 
1.2.1 主 色 描述 符 

DCD 是 MPEG-7 标准 认可 的 一 种 颜色 描述 符 。DCD 将 
图 像 的 颜色 分 为 若干 个 分 区 , 称 为 “ 粗 分 区 ”假设 同一 个 分 


区 中 的 所 有 点 均 为 相似 点 ， 分 区 中 心 是 该 分 区 中 所 有 像素 颜 
色 的 平均 值 ， 具 体 计算 为 下 式 : 

C= pa? 

~ (1) 


其 中 :Pi 表示 第 i 个 分 区 。 本 文 在 RGB 颜色 域 提取 DCD 特征 ， 
假设 每 个 像素 均 属于 一 个 分 区 ， 然 后 将 该 像素 的 颜色 值 蔡 换 
为 其 分 区 的 中 心 值 , 由 此 量化 了 图 像 的 颜色 。 对 于 每 个 像素 ， 
计算 像素 窗口 的 DCD 值 ,将 这 些 DCD 值 组 成 特征 向 量 。 图 
1 所 示 是 原 图 像 及 其 RGB 颜色 空间 中 的 DCD 图 像 。 

1.2.2 局 部 颜色 直方 图 的 特征 

直方 图 是 图 像 的 离散 统计 概率 密度 。 遥 感 图 像 一 般 为 
RGB 彩色 图 像 格 式 ， 每 个 颜色 通道 的 范围 为 [0,255]。 如 果 考 
虑 每 个 直方 图 的 全 部 256 个 bn， 那 么 计算 复杂 度 较 高 。 本 
文 将 每 个 通道 的 颜色 bin 量化 为 8 个 相等 的 bin, 考虑 每 个 像 
素 的 窗口 ， 计 算 窗口 的 局 部 直方 图 。 最 终 ， 图 像 的 每 个 像素 
被 映射 为 一 个 24 维 的 向 量 。 

1.2.3 颜色 统计 的 特征 
每 个 像素 窗口 的 第 


、 第 二 矩阵 分 别 表示 了 平均 值 与 标 
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准 偏差 。 将 像素 每 个 RGB 通道 的 平均 值 与 标准 偏差 考虑 为 
颜色 的 统计 特征 。 


(a) 原 遥感 图 像 


(b) 对 应 的 DCD 图 像 
图 1 原 图 像 及 其 RGB 颜色 空间 的 DCD 图 像 
Fig.1 Original image and DCD image in RGB color space 


2 “分 类 器 的 选择 


2.1 支持 向 量 机 (SVM) 分 类 器 

遥感 图 像 存 在 较 多 的 像素 类 别 。 将 提取 的 特征 输入 分 类 
器 来 训练 模型 的 参数 ， 对 于 每 个 像素 类 别 ， 训 练 一 个 一 对 多 
的 二 值 SVM( 支 持 向 量 机 )。 假 设 主 分 类 的 像素 为 正 样本 ， 其 
他 所 有 分 类 的 像素 为 负 样 本 ， 主 分 类 的 像素 数量 远 低 于 其 他 
分 类 的 像素 数量 ， 因 此 会 导致 SVM 模型 不 平衡 ， 并 且 偏 向 
负 样 本 。 为 了 解决 该 问题 ， 从 负 样 本 集合 中 随机 选择 个 样 
本 作为 SVM 的 负 样 本 训练 集 ， 从 而 实现 正 负 样本 的 平衡 。 
图 2 所 示 是 SVM 分 类 器 训练 数据 平衡 后 的 数据 实例 。 


一 


图 2 正 负 样本 平衡 的 SVM 训练 方法 


Fig.2 SVM training method for positive and negative sample balance 
2.2 极限 学 习 机 (ELM) 分 类 器 

ELM 原本 是 针对 单 隐 层 前 向 反馈 神经 网 络 (SLFN) 而 提 
出 的 分 类 器 ， 对 于 广义 的 SLFN( 单 层 前 馈 神 经 网 络 )，ELM 
的 输出 为 


f= Bh O) 


其 中 :pi 为 第 i 个 隐 层 节点 与 输出 节点 之 间 的 权重 , hi(x) 为 第 i 
个 隐 层 节点 的 输出 。 隐 层 节 点 的 输出 函数 hi 加) 一般 是 一 个 非 
线性 的 分 段 连续 函数 ， 例 如 : Sigmoid 函数 、Hardlimit 函数 、 
高 斯 函数 与 多 元 函数 。 

前 向 反馈 神经 网 络 的 目标 是 实现 较 小 的 训练 误差 ， 因 此 
ELM 同时 考虑 了 较 小 的 训练 误差 与 较 小 的 输出 权重 范 数 。 权 
重 范 数 越 小 ,模型 的 泛 化 性 能 越 好 。ELM 也 可 以 表示 为 以 下 
的 优化 问题 形式 : 


min:|Bl, + CIHB -TI (3) 


其 中 :ol 与 02 均 为 正 数 , p, gq=0,0.5,1.2,... 瑟 表示 有 个 方程 的 矩 
阵 ， 称 为 隐 层 输出 矩阵 ，7 为 训练 数据 的 目标 和 矩阵。 本 文采 
用 ELM 对 像素 特征 向 量 集 进 行 训练 ， 然 后 对 其 他 像素 的 特 
征 向 量 进行 分 类 处 理 。 本 文 使 用 MATLAB 中 的 ELM 工具 箱 
实现 ELM 分 类 器 。 


3 ”特性 选择 
特征 选择 是 一 个 从 个 特征 中 选 出 m 个 特征 的 离散 优化 


题 。 使 用 下 式 计 算出 所 有 可 用 子 


Ds) = +r/ Dt...+n/n) =2" (4) 


使 


中 
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的 数量 : 统 的 ACO 特征 选择 算法 中 ， 每 个 蚂蚁 基于 每 条 路 径 的 效率 


选 出 最 优 的 特征 。 假 设 特征 向 量 共 及 个 特征 ， 为 每 个 蝎 蚁 
随机 分 配 一 个 初始 化 特征 ， 第 2 个 特征 需要 训练 n-1 次 的 模 


采用 蚁 群 优化 算法 的 特征 选择 


FE 向 量 的 维度 ，s 表示 当前 特征 子 集 的 规模 。 型 训练 ， 第 3 个 特征 需要 n-2 次 的 模型 训练 。 最 差 的 情况 是 


每 次 和 迭代 、 每 个 蚂 蚊 需 完成 n(n-1)/2 次 的 模型 训练 ， 因 此 计 


以 群 砚 食 的 


: 怠 


蚁 群 优 化 算法 是 一 种 基于 迭代 与 概率 的 启发 式 算法 。 蚁 算 成 本 极 高 。 本 文 提出 一 种 新 的 方案 ， 在 保持 高 分 类 性 能 的 
自然 行为 ， 蚁 群 在 图 
E， 在 经 过 的 每 条 路 径 上 释放 信息 素 ， 根 据 信息 素 寻找 最 ”要 评估 所 有 未 访问 的 特征 ， 该 步骤 的 计算 成 本 极 高 。 从 诸多 


中 的 节点 之 间 ”同时 ， 降 低 计算 成 本 。 原 始 蚁 群 优化 算法 的 每 次 迭代 中 ， 需 


短 的 路 径 。 路 径 上 的 


言 息 素 会 随 着 时 间 而 分 解 ， 昌 蚁 则 选择 ”文献 的 实验 结果 可 看 出 ,实际 上 无 须 所 有 的 蚂蚁 执行 该 步 又。 
素 浓度 最 高 的 路 径 ， 图 中 的 最 短路 径 也 具有 较 多 的 信息 。 为 了 解决 该 问题 ， 本 文 将 特征 分 为 13 个 分 组 ， 如 表 1 所 示 。 


下 距离 的 问题 ， 图 中 日 


的 概率 可 讨 
a OF nt 


PD=4 ,nr, (OF 


以 选择 的 概率 也 更 高 。 将 特征 选择 问题 建 模 为 ” 表 中 LH_meanL 表示 小 波 LH 频带 工 分 量 的 范 数 均值 ， 


0, 其 他 | 
:天 表示 当前 尚未 被 访问 的 特征 


关联 的 信息 素 值 与 
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下 式 : 色 ，DCDR;G 表示 G 颜色 分 量 的 主 色 ,DCDRIiB 表示 B 颜色 
et 分 量 的 主 色 ，Rn&R; 量 表示 R 颜色 分 量 的 均值 与 标准 偏差 ， 

(5) ”Gmn&G; 量 表示 G 颜色 分 量 的 均值 与 标准 偏差 ,Bm&B; 量 表示 

况 B 颜色 分 量 的 均值 与 标准 偏差 ，R; 表 示 R 颜色 分 量 的 直方 图 

着 ，Ti(D 与 wy 分 别 为 特征 。 bin，Gi 表 示 G 颜色 分 量 的 直方 图 bin，Bi; 表 示 B 颜色 分 量 的 


启发 值 , y 与 5 分 别 为 信息 素 与 启发 信息 。 直方 图 bin。 


的 节点 表示 特征 ， 边 表示 选 。 LH_meanH 表示 小 波 LH 频带 了 HH 分 量 的 范 数 均值 , HL_meanH 
Et151。 蚁 群 在 图 中 游 走 之 后 ， 选 出 最 优 的 特征 。 表示 小 波 HL 频带 H 分 量 的 范 数 均值 ，HL_meanL 表示 小 波 
息 素 、 启 发 值 与 节点 (特征 ) 关 联 ， 蚁 群 大 在 ” HL 频带 LL 分 量 的 范 数 均值 DCDRiR 表示 R 颜色 分 量 的 主 
| 算 为 


时间 ft 时 曙 蚁 在 特征 i 释放 的 信息 素 计算 为 下 式 : 表 1 小 波 特征 与 颜色 特征 的 分 组 
y. ( FC- | SE 0) Table ] Grouping of wavelet features and color features 
St 分 组 编号 此 征 意义 分 组 编号 。 ”特征 意义 
Ar(D = if (ie s(t) (6) 1 LH_meanL( 小 波 特征 ) 8 Ron&R:( 颜 色 特 征 ) 
2 LH_meanH( 小 波 特征 ) 9 Gm&Gs( 颜 色 特 征 ) 
0 3 HL_meanH( 小 波 特征 ) 10 。。 轧 &&B:( 颜 色 特 征 ) 
P:FC 表示 所 有 特征 的 数量 ，S*(D) 与 1S*D)| 分 别 为 蚂蚁 在 4 HL_meanL( 小 波 特征 ) 11 Ri: 颜色 特征 ) 
征 子 集 与 特征 长 度 。 瑟 (SOD) 表 示 对 子 集 SK(7) 5 DCDR 让 (颜色 特征 ) 12 Gx( 颜 色 特 征 ) 
能 的 评估 结果 ，4Y 与 y 分 别 是 平衡 分 类 器 性 能 与 特 6 DCDRiG( 颜 色 特 征 ) 13 Bi( 颜 色 特 征 ) 
度 的 权重 ， 节 点 的 信息 素 更 新 方法 如 下 式 ; 7 DCDRiB 
ee 特征 分 组 的 优点 是 每 个 蚂蚁 仅 需要 评估 每 组 的 一 个 特 
妇 征 。 因 为 每 个 分 组 内 的 特征 对 于 分 类 准确 率 具有 相似 的 效果 ， 
量 , p 为 信息 素 的 挥发 速率 ，8 为 当前 的 最 。 所 以 无 须 对 每 个 分 组 的 所 有 特征 进行 训练 。 例 如 : 如 果 


了 最 优 蚂 蚁 对 节点 的 信息 速度 影响 大 于 其 。 LH_meanL 是 在 第 上 次 迭代 从 分 组 1 中 选择 的 特征 ， 那 么 特 
此 之 后 的 近代 中 凡 群 会 集中 于 最 优 解 。 基 。” 征 LH_stdL 将 是 第 kt1 次 迭代 的 一 个 候选 解 ， 说 明 该 分 组 中 
征 选 择 方案 如 图 3 所 示 。 的 所 有 特征 均 有 机 会 被 选择 。 
四 第 一 个 特征 选择 方案 中 ， 同 时 对 所 有 的 像素 分 类 进行 特 
4 征 选 择 处 理 ， 特 征 选择 方案 通过 最 大 化 所 有 像素 的 分 类 准确 
| 率 ， 提 取出 最 优 的 特征 。 因 此 该 方法 中 各 个 分 类 选择 相同 的 
特征 子 集 。 
Sa 第 一 个 特征 选择 方法 的 步骤 如 下 所 示 ; 
HE) a) 设 置 蚁 群 的 规模 。 将 每 个 特征 与 图 中 每 个 节点 建立 映 
射 ， 为 每 个 节点 设置 一 个 随机 的 信息 素 水 平 。 设 置 最 大 迁 代 
次 数 与 结束 条 件 。 
由 .并 b) 将 特征 集 分 组 。 
KARL 0) 每 个 蚂蚁 创建 一 个 解 ， 每 个 蚂蚁 从 初始 解 (初始 节点 ) 
R 信息 素 最 高 的 路 径 。。 羽 &- 开始 游 走 ， 并 且 将 该 解 标记 为 “已 访问 ”。 
图 3 基于 ACO 的 特征 选择 方案 d) 每 个 蚂蚁 从 每 个 特征 分 组 随机 选择 一 个 特征 ， 并 将 该 
Fig. 3 Feature selection scheme based on ACO 特征 保存 于 队列 a。 使 用 式 (5) 计 算 a 中 所 有 特征 的 概率 值 ， 
3.2 本文 的 特征 选择 方案 概率 值 最 高 的 特征 标记 为 “已 访问 ”。 
和 相关 性 最 高 的 特征 ， 不 仅 能 够 降低 特征 向 量 的 维 e@ 如 果 蚂 蚁 未 能 满足 阔 值 条 件 ( 式 (8)， 那么 返回 步 又 d)。 
度 与 分 类 器 的 计算 成 本 ， 而 且 能 够 提高 系统 的 分 类 准确 率 。 
本 文 的 像素 分 类 系统 中 ， 在 特征 提取 阶段 完成 特征 选择 的 处 A =p*exp 上 Orexp (8) 
四 其 中 :EN 是 蚁 群 目前 已 选择 特征 的 特征 基数 ，N 是 特征 的 总 
所 有 分 类 使 用 同一 特征 集 (第 一 种 特征 选择 方案 ) 数量 ,9 与 @ 分 别 是 控制 特征 规模 与 F_measure 效果 的 参数 ， 
甘于 ACO 的 特征 选择 方案 是 一 种 基于 封装 的 方案 , 传 gt+w=1。 最 终 ， 所 有 蚂蚁 成 功 获得 相应 的 特征 子 集 。 
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考虑 了 三 条 信息 素 的 更 新 规则 : (a) 每 个 蚂蚁 在 特征 上 释 
放 一 些 信息 素 ; (b) 一 些 信息 素 随 着 时 间 而 挥发 ，(o) 最 优 的 照 
蚁 对 信息 素 更 新 具有 决定 性 效果 ， 从 而 其 他 蚂蚁 可 以 跟随 优 
质 蚂蚁 的 路 径 。 
人 ) 蚁 群 的 所 有 蚂蚁 执行 步 又 c)~e)。 
多 每 个 蚂蚁 在 其 路 径 中 的 特征 上 释放 信息 素 ， 并 将 该 节 
点 标记 为 “已 访问 ” 信息 素 更 新 方法 计算 为 下 式 ， 
Apher(i,k) = .( Acc_model(k)) 
+ 人 [于 cc (9) 


F_C 


其 中 :Apherl9 是 曲 蚁 上 在 特征 i 上 释放 的 信息 素 ， 
Acc_model0 表 示 蚂 蚁 提取 特征 所 训练 的 模型 分 类 准确 率 。 
F_N 与 已 MD 分 别 是 特征 总 数量 与 蚂蚁 K 路 径 中 的 特征 数 
量 。 最 终 ，v 与 8 分 别 控制 分 类 准确 率 与 特征 子 集 长 度 ， 并 
且 x+p=1。 

hb) 搜索 F-measure 值 最 高 的 蚂蚁 。 

j 更 新 信息 素 : 


pherlt +1)= (1- p).pher(t) 


+》 Apher(i,k) + Apher(i, Antina) 
k=] 


(10) 


删除 己 有 的 蚂蚁 ， 随 机 产生 新 蚂蚁 。 
区 如果 满足 结束 条 件 ， 则 执行 步骤 D， 和 否则 跳 至 步骤 c)。 
DF-measure 值 最 大 的 路 径 即 为 最 优 解 。 

图 4 所 示 是 第 一 种 特征 选择 方案 的 流程 框图 。 


全 部 特征 集 四 
了 一 达到 改 群 阅 值 ? = 


订 妈 化 和 人、 节点 信 息 素 、 结 | Bs 
一 所 以 如 坟 结 束 游 老 ? 一 一 


特征 分 组 是 
T 计算 所 有 蚂蚁 的 delta 信 息 素 
建立 蚁 群 + y 
T 搜索 分 类 性 能 最 优 的 蚂蚁 
选择 蚁 群 的 初始 特征 - 了 
更 新 信息 素 
从 每 个 种 群 随机 选择 一 个 特征 ， 并 | | 
且 保 存 于 R 变 量 = 达到 结束 条 件 ? 一 


-证 
保存 分 类 性 能 最 优 的 特征 集 


Y 
基于 概率 从 R 中 选择 一 个 特征 


图 4 第 一 种 特征 选择 方案 的 流程 
Fig.4 Flow of the first feature selection Scheme 
3.2.2 不 同 分 类 使 用 不 同 特征 集 (第 二 种 特征 选择 方案 ) 
第 一 种 特征 选择 方案 将 选择 的 最 优 特征 集 作 为 7 个 像素 
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类 的 最 优 特征 子 集 。 第 一 个 方法 的 特征 降 维 阶段 ， 研 究 了 像 
素 分 类 的 主要 属性 ， 并 且 提 出 了 包含 小 波 特 征 与 颜色 特征 的 
特征 向 量 。 对 7 个 分 类 分 别 使 用 不 同 的 最 优 特 征集 ， 可 提高 
分 类 的 准确 率 ， 图 5 所 示 分 别 是 像素 分 类 7 与 像素 分 类 4 的 
图 像 ， 图 中 小 波 特 征 的 水 平 、 垂 直 频 率 是 分 类 7 的 不 相关 特 
征 ， 而 是 分 类 4 的 高 度 相关 特征 。 


(a) 像 素 分 类 7 的 图 像 实例 


(b) 像 素 分 类 4 的 图 像 实例 
图 5 像素 分 类 7 与 像素 分 类 4 的 图 像 
Fig.5 Pixel classification 7 and pixel classification 4 
第 二 种 特征 选择 方法 为 每 个 像素 分 类 分 别提 取 了 一 个 最 
优 特 征 子 集 ， 虽 然 该 方法 的 计算 时 间 较 长 ， 但 是 该 程序 为 线 
下 程序 ， 可 以 保证 线 上 程序 的 实时 性 。 图 6 所 示 是 第 二 个 特 
征 选择 方案 的 流程 框图 。 
这 种 类 型 的 特征 选择 方案 为 每 个 每 个 像素 分 类 设置 一 个 
最 优 的 特征 子 集 ， 从 而 保证 每 个 像素 分 类 有 具有 最 相关 的 特征 
子 集 ， 以 期 提高 总 体 的 像素 分 类 准确 率 。 第 二 种 特征 选择 方 
案 的 步骤 如 下 所 示 : 
a) 每 个 像素 类 运行 步骤 b)~d)。 
b) 运 行 第 一 种 特征 选择 方案 。 
c) 分 别 计算 每 个 分 类 的 最 优 特征 子 集 。 
d) 分 别 训练 每 个 像素 分 类 的 分 类 器 。 


全 部 特征 集 


基于 蚁 群 优化 群 优化 


的 特征 选择 Ey 蚁 群 
采用 ELM 的 


王 
Nr 


像素 分 类 器 


特征 子 集 | 人 下 
的 


局 更 新 信息 素 


第 一 个 类 的 最 优 特征 
子 集 


子 集 


第 二 个 类 的 最 优 特征 


蚁 群 
ed 2 HJ jy 畦 征 子 集 | “了 咯 
i | 好 
更 新 信息 素 
第 七 个 类 的 最 优 特征 
子 集 


图 6 第 二 个 特征 选择 方案 的 流程 框图 


Fig.6 Flow ofthe second feature selection Scheme 
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4 ”实验 环境 与 参数 设置 


采用 Mars Exploration Rover Spirit 数据 集 作 为 
benchmark 数据 集 ， 该 数据 集 共有 94257 个 图 像 ， 图 像 大 小 
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均 为 512x512。 图 7 所 示 是 该 数据 集 的 三 个 图 像 实 例 。 文 献 


0 将 该 数据 集 分 为 7 个 分 类 ， 如 图 8 所 示 ， 将 文献 上 的 分 类 
结果 作为 先 验 知识 。 


(a) 


图 7 遥感 


(b) 


(c) 


图 像 的 实例 


Fig.7 Examples of remote sensing images 


(a) 第 一 类 (岩石 轨迹 ) 


人 ) 第 二 类 (小 的 黑色 石头 与 沙 ) 


(e) 第 五 类 (波浪 岩石 ) 


图 8 


本 算法 的 编程 环境 为 : Intel CPU， 主 频 为 2.26 GHz， 内 
存 为 6 GB。 每 个 像素 的 窗口 规模 为 21x21， 计 算 灰 度 图 像 每 
个 窗口 的 小 波 分 解 ， 并 且 计 算 窗 口中 行 、 列 的 范 数 ， 最 终 每 
个 像素 的 特征 向 量 均值 与 标准 偏差 共有 8 个 小 波 特征 。 颜 色 
特征 的 主 颜 色 量 设 为 8， 因此，R、G、B 三 色 通 道 共有 24 
个 特征 ， R、G、B 通道 的 均值 与 标准 偏差 组 成 6 个 颜色 的 统 
计 特 征 。 使 用 R、G、B 三 个 直方 图 的 4 个 bin， 共 计算 出 12 
个 特征 。 最 终 ， 每 个 像素 均 有 一 个 50 维 的 特征 向 量 。 

本 实验 的 7 个 分 类 共有 7 个 ELM 分 类 器 ， 每 个 分 类 对 
应 一 种 地 表 类 型 。 在 每 个 分 类 器 的 训练 过 程 中 ， 将 对 应 的 地 
表 类 型 作为 正 样本 ， 其 他 像素 作为 负 样 本 。 表 2 所 示 是 实验 
中 所 有 7 个 分 类 的 正 、 负 样本 数量 ， 每 组 数据 集 的 60% 像 素 
作为 训练 集 ， 剩 下 的 40% 作 为 测试 集 。 
表 2 Mars Exploration Rover Spirit 数据 集 7 个 分 类 的 正 负 样本 数量 

Table 2 7 classifications in mars exploration rover Spirit data set 

分 类 1 2 3 4 5 6 7 


正 样 本 5651 5470 4941 10771 8682 4959 5475 
负 样 本 5734 4365 5412 12544 9322 4566 7165 
样本 总 量 11385 9835 10353 23515 18004 9525 11640 


KNN、GA 与 本 算法 在 不 同 的 参数 设置 下 分 别 进行 了 实 
验 ， 测 试 不 同 参 数 下 三 个 算法 特征 选择 与 像素 分 类 的 性 能 。 
分 类 器 模型 采用 SMO 算法 (sequential minimal optimization) 
作为 SVM 的 学 习 方法 , 最 大 办 代 次 数 为 20000， 内核 缓存 约 
束 为 1000。 多 层 感 知 机 (multi layer perceptron) 的 范围 为 
[-0.01,0.01]， 高 斯 核 为 线性 核 函数 ， 并 且 o=4。 


(上 第 六 类 (大 岩石 ) 
遥感 图 像 的 7 个 分 类 实例 


Fig.8 classification examples of remote sensing images 


(c) 第 三 类 (中 等 石头 与 沙 ) (d) 第 四 类 (分 层 岩石 ) 


(g) 第 七 类 (平坦 岩石 ) 


采用 几 个 常用 的 分 类 器 性 能 指标 评估 本 算法 的 性 能 ， 分 
别 为 精度、 召回 率 、F-measure、 准 确 率 ，F-measure 定义 
为 


F-measure =2x(Precision x Recall)Precision + Recall 


其 中 precision 为 精度 ，recall 表示 召回 率 。 
5 ”实验 结果 与 分 析 


第 一 个 特征 选择 方案 的 实验 结果 

为 了 比较 不 同 分 类 器 的 效果 ， 对 ELM，SVM 与 KNN 3 
个 分 类 器 分 别 进行 实验 分 析 ，ELM，SVM 与 KNN 分 别 采用 
13、17、18 个 特征 。 图 9 所 示 三 个 分 类 器 获得 的 分 类 性 能 ， 
图 中 可 看 出 ，ELM 分 类 器 实现 了 最 高 的 性 能 。 


(13) 


5.1 


100 
90 上 
及 80 上 
上 
怠 
3 70 上 
L —s— ELM_RBF 
eol —e— SVM_MLP 
一 KNN(k=5) 
50 LL 上 1 1 1 L 1 1 
0 2 4 6 8 10 12 14 16 18 
特征 数量 (个 ) 
图 9 三 个 分 类 器 的 分 类 性 能 (第 一 个 特征 选择 方案 ) 
Fig.9 Classification performance of 3 classifiers (1st) 
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5.2 第 二 个 特征 选择 方案 
第 二 个 特征 


等 : 基于 蚁 群 优化 与 独立 特征 集 的 路 感 图 像 实 


选择 算法 为 每 个 分 类 器 提取 了 一 个 优化 的 特 
征 子 集 。 对 ELM 于 SVM 两 个 分 类 器 分 别 进行 了 实验 ,结果 


如 图 10 所 示 。 从 结果 可 看 出 , 极限 学 习 机 的 特征 数量 略 多 于 
SVM 分 类 器 。 
24 上 | —=— ELM_ RBF 
—e— ELM_line 
22 上 | 一 一 SVM_MLP 
—v— SVM _line 
20 上 
Cel 
晒 18 
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图 11 所 示 是 ELM 与 SVM 分 类 准确 率 的 结果 ， 图 中 可 
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Fig. 12 Computing time of remote sensing image classification 
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Fig. 13 Examples of image segmentation results 
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